Towards a picosecond transform-limited 
nitrogen- vacancy based single photon 

source 



Chun-Hsu Su, Andrew D. Greentree, and Lloyd C. L. Hollenberg 

Quantum Communications Victoria, School of Physics, 
The University of Melbourne, Victoria 3010, Australia 
chsu @ ph. unimelb. edu. au 



Abstract: We analyze a nitrogen-vacancy (NV^) colour centre based sin- 
gle photon source based on cavity Purcell enhancement of the zero phonon 
line and suppression of other transitions. Optimal performance conditions of 
the cavity-centre system are analyzed using Master equation and quantum 
trajectory methods. By coupling the centre strongly to a high-finesse opti- 
cal cavity [Q - ~ 10^), V - A^] and using sub-picosecond optical 
excitation the system has striking performance, including effective lifetime 
of 70 ps, linewidth of 0.01 nm, near unit single photon emission probability 
and small [^(10^^)] multi-photon probability. 

OCIS codes: (230.0230) Optical devices; (230.6080) Sources; (160.4760) Optical properties; 
(160.2220) Defect-center materials 
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1. Introduction 

A source that produces single photons on demand is an invaluable tool for precision optical 
measurement [1,2] and is a crucial building block for many quantum computing and com- 
munication applications. For an attenuated laser, the number of photons per pulse follows a 
Poisson distribution and the multi-photon probability becomes negligible only in the limit of 
small mean photon number at the expense of the single-photon probability. In quantum comput- 
ing, using single photons to store and transport quantum information is natural as information 
can be easily encoded and manipulated over photonic degrees of freedom e.g. polarization. In 
linear optical quantum computing (LOQC) [3,4], it is straightforward to perform single qubit 
operations on photons with elementary optical components and projective measurements gen- 
erate photon-photon interactions. In quantum communication, single photon sources can be 
used for unconditionally secure quantum key distribution (QKD) protocols [5,6]. To date, sin- 
gle photon generation has been demonstrated with a variety of single quantum emitters such as 



atoms [7-10], ions [11], molecules [12, 13], diamond colour centres [14-17] and semiconductor 
quantum dots [ 1 8-2 1 ] . 

Diamond defects hold promise as a platform for solid-state quantum optical and quan- 
tum computing applications [22], and for the study of condensed-matter analogues [23-25]. 
One of the most well-studied systems, garnering much attention lately, is the optically active 
negatively-charged nitrogen-vacancy defect (NV^ centre). The centre consists of a substitu- 
tional nitrogen atom and an adjacent vacancy in the carbon lattice. It forms naturally or may be 
engineered within a diamond matrix using techniques such as single ion implantation [26-28] 
or chemical vapour deposition [29]. It has a combination of remarkable properties that render 
it a suitable single photon source candidate. These include robustness against photobleaching, 
structural stability at room temperature and demonstrated antibunching, which is the hallmark 
of a single photon source [14]. The centre has also been used to realize Wheeler's delayed- 
choice experiment [30]. However, the centre has a relatively long photoluminescence lifetime 
of 1 1.6 ns and broad spectral width of 150 nm [31], which is not optimal for daylight or opti- 
cal fibre operation in QKD [32, 33]. Furthermore, the photons are not time-bandwidth limited 
(or indistinguishable) for the purpose of LOQC where photon indistinguishability is crucial for 
Hong-Ou-Mandel interference [34] and hence quantum gates [35]. 

Preparing the centre in a high-finesse quantum cavity offers a solution to these problems. 
Cavity quantum electrodynamics has been shown to induce a single-photon Kerr nonlinear- 
ity [36] and assist quantum gate operation [37,38]. Under strong photonic confinement, the 
quantum emitter (the centre in this case) in the cavity interacts coherently with photon states 
with the effect of modifying photon-emission dynamics. As a result, we show that the spectral 
properties of the centre can be improved to fulfill the stringent criteria for quantum information 
applications. Additionally, the emission can be directed into an application or experiment as 
desired. A suitable cavity is the planar photonic-band-gap (PBG) cavity that defines an excel- 
lent cavity with small mode volume (of order one cubic wavelength) and low loss that provides 
strong centre-photon coupling [39^1] suitable for our purposes. 

Advances in fabrication techniques are nearing the stage where preparing a single crystal 
diamond with PBG structures using lithography and lift-off [42^6] and placing an individual 
NV centre in the centre of the cavity using ion implantation techniques may be possible. The 
latter technique permits locating the centre to achieve full emission enhancement. We note that 
cavities have been used to enhance the emission of quantum-dots [18,20] and atoms [10] for 
photon generation. Our study is made in a similar spirit, but a diamond-based device has the 
advantages of robustness against the environmental noise combined with simplicity of its setup. 
Here, without loss of generality, we theoretically study the effects of placing the centre within 
a high-finesse single-moded diamond PBG cavity. We establish the operation criteria (cavity 
specification and excitation scheme) for an efficient cavity-centre based single photon source. 

2. Theory 

A model of an NV^ centre has been proposed as a vibronic system with ground (g) and excited 
(e) electronic states, each given by a series of vibrational sublevels \gi) and \ej) respectively, 
where / and j label the vibrational states [47]. Its emission spectra consists of several phonon 
lines /PL, corresponding to the transition |eo) — \gi)- For simplicity, |eo) is denoted as |e) from 
now on. Phononic relaxation of the excited state is much faster than the radiative relaxation 
to the ground state and the photonic transition probabilities were calculated by Davies and 
Hamer [47] under the WKB approximation. Because of the rapid phononic transitions, we 
treat the centre as an atomic system with a single excited state \e) and a ground state with ten 
sublevels {|^,)} (Fig. [Til. The fluorescence from the centre corresponds to a transition from the 
excited spin triplet state {^E) to an electron spin triplet ground state (^A) and the dynamics are 



Fig. 1. Theoretical model of a cavity-centre system for single photon generation: the NV 
centre is modelled as a multi-level atom with a single excited state \e) and a ground state 
with vibrational sublevels {|g))}. The centre is pumped with an external classical field 
r{t) (white arrow) acting as the trigger pulse, the transition \e) — |g,) is coupled to a lossy 
single-modal cavity with coupling strength (grey) and cavity decay rate K (black), jg.^ 
(dashed) is the atomic decay rate for radiative transitions \e) — \gi) while yg,„i;„ (dotted) are 
that for the non-radiative phononic transitions \gn) — Igm)- 

influenced by the presence of a possible metastable state ('A). We have explicitly ignored the 
metastable state as its role in de-excitation process remains unclear [31,48,49]. Also, in the 
absence of strain and magnetic field, the transitions from the = spin sublevel of the ground 
state are found to be spin-conserving [31], hence we assume this in our treatment. 

To study the transient behaviour of the cavity-centre system, we extend the basic Jaynes- 
Cummings (JC) model [50,51], which treats the interaction between a single-mode electro- 
magnetic field or cavity of resonance frequency tOc and a two-level atom, to consider our more 
complex atomic system. In the rotating-wave approximation (RWA), the JC Hamiltonian of our 
cavity-centre system is expressed in terms of the atomic projection operators dap = | o:) (j3 1 and 
the annihilation (creation) operators d {a') for single cavity mode as, (h = I) 

Na-2 ^ Na-2 

jr-"^ = ^ (Og. + (Oedee + OcO^ d + X ^ (^''^^ + h-C) , (1) 

(=0 ^ 1=0 

where A'a = 11 is the total number of atomic states, Oa is the energy of the atomic level \a) 
and is the cavity-centre coupling constant between atomic transition \e) — \gi). In the dipole 
approximation, the coupling is i2, = di[(Oc / {'2-TieQV)Y 1^ where di is the dipole moment of the 
respective transition. The evolution of the cavity-centre system obeys the Liouville equation of 
motion for the density matrix p, 

do ^'^'^ ^''^^ 

i=o i=0 

with the Lindbladian terms for some operator O, 

^[d,p] = dpd'' ~^{d''dp+pd''d) . (3) 

The spontaneous transition (|e) — \gj)) couples to any non-cavity field modes at the character- 
istic rates Ygje and non-radiative phononic decays from \gi+i) to \gi) with rates Ygigj^^ - Here 
b'^ represents the creation operator for electromagnetic (waveguide) mode outside the cavity. 



which the cavity couples to via decay rate K. The rate K = C0c/{2Q) is parameterized by the 
quality factor of the cavity Q. Incoherent excitation, acting as the trigger for photon emission, 
is represented by the phenomenological term with a pump absorption rate r(f ). In this model, 
we have explicitly ignored thermal broadening by assuming a zero temperature operating en- 
vironment. Broadening can be introduced phenomenologically for a more realistic estimate of 
the linewidth of the emitted wave packet. 

Efficient single photon generation requires minimizing loss and fast outcoupling of the exci- 
tation via the cavity channel. The appropriate regime to optimise the output is the strong Purcell 
regime [52], K > Q.i^ y^.g, where the rate of coherent coupling between the centre and cavity 
mode, n?/fc, dominates the rates j^.^ (Vy) of the incoherent coupling to the non-cavity modes. 
The cavity mode is chosen to match the transition \e) — The cavity loss rate fc sets the time 
scale for photon outcoupling and when greater than H,, suppresses the vacuum Rabi oscilla- 
tions, which otherwise lead to unwanted spectral features on the output photon. For a two-level 
atom in a cavity, altogether embedded in a medium of refractive index «, the enhancement of 
emission into the cavity is parameterized by the Purcell factor, the ratio of the emission rate to 
the cavity mode to the unmodified rate into the non-cavity modes [53-55], 



3 (Ac/n)^ 6 
4;r2 V" 



Fp= (4) 



where V is the cavity mode volume and Ac is the cavity wavelength. 

To determine the fraction of pulse cycles that lead to useful output, we must consider the 
fraction of emitted photons from the atom that enter the desired cavity mode. This fraction 
is set by the spontaneous coupling factor j3 = Fp/{1 +Fp). We demand j3 to be near unity, 
which implies a high-Q/V cavity to maximize efficiency, but a small Q for fast cavity loss. 
These considerations lead to an upper bound for Q, or equivalently, a lower bound on K given 
by fc > 2Q.i. Beyond this, we enter the strong cavity regime [54]. In this regime, the vacuum 
Rabi oscillations are not sufficiently suppressed by cavity decay leading to an effective timing 
jitter or equivalently temporal/spectral features which will degrade overall device performance. 
Alternatively, if fc < y^ije for large Q, the excitation wiU be outcoupled as atomic decoherence, 
which leads to loss of photons through non-guided modes, and manifests as an increased zero 
photon probability. A rigorous treatment of cavity-assisted emission was made in Ref. [56]. 



3. Analysis and results 

To investigate dynamical processes of photon generation from the cavity-centre system de- 
scribed by Eq. |2] we use direct numerical integration to study the specifications on the trig- 
ger pulse that ensures high-fidelity single photon emission and the effect of the cavity on its 
characteristics. Additionally, we use the quantum trajectory approach to simulate photodetec- 
tion experiments [57,58]. We show that it yields results that agree with the former approach 
and demonstrate bit-stream photon generation and, most importantly, antibunching with the 
Hanbury-Brown-Twiss (HBT) experiment. 

3.1. Determination of the pulsed excitation parameters 

The specification of the pulsed excitation scheme for the cavity-centre system is crucial to 
ensure efficient and high-fidelity single photon emission. The complete treatment of Eq. |2] 
for this study is computationally demanding as it requires adopting a considerably large state 
space given by {|a)}atomic ® {|0c), | Ic), |A'c)}cavity ® {|0w), |lw), lA^w)} waveguide, where 
the state evolution is taken to span over A^w (or equivalently Nc) excitation states. 'C denotes 
the field mode in the cavity and W labels the external waveguide mode. However, for the 
purposes of understanding just the photon emission, we may reduce the system state space by 
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Fig. 2. a. Probability of the cavity-centre system (coc = COzpl^ V = ^zpO '■° ^ single 
photon per top-hat excitation pulse as a function of pulse width T and absorption rate tq. 
Dash-dotted line denotes the pulse width parameter used in Ref . [59] where the illumination 
irreversibly transforms the centre into a different centre and is therefore a practical cutoff. 
Circle labels the parameters [yield Pi = 0.996 and P>2 ~ i^(10^^)] used for the demon- 
stration of single-photon generation with the cavity-centre system, b. Zero (dotted), single 
(solid) and multi- (dashed) photon probability as a function of pulse width with constant 
absorption rate ro = lO'"' Hz. 

setting Nc = 2, Nw = and enforcing a choice of parameters for the pulsed excitation that 
limits system population to the < 2-quantum states throughout the cycle. 

For the excitation, we assume a top-hat form of the trigger pulse r{t) = rg for t e [0,r], 
where ro is the absorption rate and T the pulse width which must be much shorter than the 
photoluminescence lifetime. In the strong Purcell regime the centre may be approximated as a 
two-level system in resonance with the cavity, and for very short times we ignore the effects of 
spontaneous emission and cavity outcoupling. Under these approximations the lower bound for 
zero photon emission probability (Pq) and the upper bound for one (Pi) are 



and upper bound for multi-photon probabilities is P>2 = 1 ^ (^b + ^'i), where rj = {rQ — 
16n?)'/^. In the limit ro ^ ii,-, Eq. |6]reduces to Pi = 1 — exp(— roT), but the second order 
term reveals that the single photon probability decreases with increasing pulse width accord- 
ing to exp(— 4n?r/ro) — Pq, where Pq = [2 — exp(— 4£2?r/ro)]exp(— roT). There is therefore 
a clear trade-off between the requirement to produce a photon on demand and that of having 
no more than one photon per pulse. Fig. |2]illustrates the single photon emission probability for 
an NV centre in a cavity resonant with the zero-phonon line (ZPL) transition (|e) — \go)), as a 
function of excitation parameters. In contrast to an attenuated laser where small P>2 is achieved 
in expense of Pi, the cavity-centre system offers Pi ~ while keeping P>2 < ff{10^^). The 
mean photon number per pulse is « = ^, /P, w Pi + 2P>2 if P, is negligible for ; > 3, for these 
parameter choices. 

The photon source must operate in the limit of large ro and maintain short T to ensure near 
unit single photon probability and efficient operation. However, Dumeige et al. showed that 
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Fig. 3. Evolution of a cavity-centre system (coc = cozpl^ V = ^^pL' ^ ~ 2-5flo) in re- 
sponse to a top-hat excitation pulse (rg = lO'^ Hz, T = 0.56 ps). a. Population in |e,0c,0iy) 
(the excited centre, black dash-dotted) and \go,Oc^ Iw) or pww (the outcoupled waveguide 
mode, black/red solid) as a function of time. b. Time derivatives of pww^ proportional to 
the output intensity, with an integrated area of 0.99 (solid) and of l.OI (dashed red). Simu- 
lation is performed by direct integration of Eq.|2](black solid/dash-dotted) and by quantum 
trajectory approach as a direct photodetection experiment (red dashed). 

under intense femtosecond illumination, the centre becomes photo-ionized, resulting in blink- 
ing [59], thus this sets the lower bound for the pulse width. 

The optimal excitation parameters for the operation of cavity-centre system depend on the 
desired application for the photon source. The source can be tailored by varying pulse parame- 
ters T and ro to optimize properties Pi and n for a specific application. We adopted technically 
feasible values T = 0.56 ps and ro = 10^^ Hz that yield Pi = 0.996 and P>2 ^ ^(10^^). 

3.2. Simulation of single photon generation with the cavity-centre system 

The analysis in Sec. 3.1 allows a determination of the parameter range for single photon emis- 
sion, and this is obviously the regime in which we want our device to work. In this section, 
therefore, we work in the single photon emission regime and assume single photon output. This 
affords a considerable saving in computational complexity and we thus truncate the photonic 
state space to one excitation by setting Nc ~ Nw = 1- We may then solve Eq.|2]numerically to 
simulate the response of the cavity-centre system to a pulsed excitation, with the result given in 
Fig. [3] The cavity of volume V ~ ^zpl chosen to maximize coupling and is in resonance with 
the ZPL. We choose Q of 36500 so that K = 2.5i2o- A photon is being issued at a mean time 
70 ps from the pump with excitation outcoupling into the external waveguide mode illustrated 
by increasing population in \go,Oc, Iw) (or pww = (go,Oc, lw\p\8o,Oc, Iw))- The integral of 
the derivative pww is near unity at 0.99, as required of a single photon pulse. Fourier transform 
of the temporal profile yields an emission spectrum centered at Xzpl with effective linewidth of 
0.01 nm. 

Due to the difference between cavity-centre coupling iio and cavity outcoupling fc, the resul- 
tant photon pulse does not take the form of a Gaussian function which is optimal for LOQC [35] . 
In principle. Stark tuning [60] can be used to optimize the atom-cavity coupling to reshape the 
temporal pulse profile [61,62] and suppress timing jitter Note that there is a gentle hump at 
300 ps, representing the Rabi remnant, that can be eliminated with greater cavity damping. 

To simulate photodetection experiments via the quantum trajectory approach, we again adopt 




Fig. 4. Comparison of the probability of the cavity-centre system ((Bf = (Ozpl, V = ^zpD '° 
emit a photon of COzpl via the cavity channel and to emit ;PL photon via atom decoherence 
as a function of cavity quality factor Q. 



a two-level model of the centre, but truncate At = 4 to allow for the possibility of multi-photon 
occupation and subsequent emission. There is a good agreement between this result and that 
from direct integration, as shown in Fig. [3] 

It is instructive to observe how the probability of the system to emit a photon via the cavity 
and via atomic decoherence varies under the influence of the cavity. We recalculate the emission 
(Fig.lH, but this time treat (A'^-l) vibrational ground states, and only consider one photonic 
excitation. In the weak Purcell regime, the excitation is outcoupled via the atomic decoherence 
channel. The relative transition probabilities of the respective phonon lines jPL approach the 
unmodified atomic branching ratios Ygje in the limit of small Q. However, in the strong Pur- 
cell regime with Q ^ ff{10'^ — 10^), the ZPL transition is enhanced by the cavity while other 
transitions are suppressed accordingly. The excitation is predominantly outcoupled via the cav- 
ity relaxation channel, representing the optimal regime for the single-photon source. We note 
that such 2's are technically achievable and have been demonstrated in Ref. [39,40] in silicon, 
and are feasible in diamond [43]. Finally, in the strong cavity regime, the Purcell enhancement 
diminishes as the time scale for cavity relaxation becomes much longer with K < y^.g. 

We have also considered the possibility of enhancing the higher order phonon side bands. 
Such an enhancement is attractive from the point of view of shifting the emission further into 
the IR, and taking advantage of the increased dipole moments. The presence of the vibrational 
sublevels and phononic decays introduces an additional complication to the dynamics of the 
basic JC model. By explicitly considering a reduced centre model with \e) and (/ = 0, 1,2 
only) that couples the cavity that is in resonance with the IPL, we have obtained an analytical 
expression to the modified decay rate by solving the master equation for its eigenfrequencies 
and assuming i2o ^ ('^+ Tsmfti)- Single excitation relaxes at an overall damping rate, 

7..,.,,^z,..,, V/(2a.r:,..-2r.., "'^'"°'^ ™ 

Eq. |4]is recovered by setting Jf^^j.^ ~ lumHn = 0. The phononic decay reduces the coherence 
between the states Xq) and |e,Oc), inhibiting the Purcell enhancement. Hence, to enhance 
the higher phonon transitions, one requires i2o of order i^(7g,„i;„) or a cavity whose modal 
volume is much smaller than the scale of wavelength^. Although technically demanding, we 
note that sub-wavelength confinement may be possible with plasmonic cavities [63]. 




Fig. 5. Photon correlation liistogram of emission from tfie cavity-centre system under 
pulsed excitation of top-hat functional form obtained using a HBT setup with quantum tra- 
jectory approach. The simulations involve a. excitation pulse of temporal width T = 0.56 ps 
and constant absorption rate vq = lO'-' Hz at a repetition rate of 1 GHz for a trajectory time 
of 5 /xs, and b. excitation pulse of varying temporal widths and constant ro = lO'"' Hz at 
repetition rate of 0.5 GHz for total time 1.5 ms. 



3.3. Antibunching with a HBT setup 

Photon correlation obtained from a HBT experiment is a test for single-photon emission. The 
HBT setup is simulated with a pulsed excitation of rate 1 GHz for a simulation time of 5 /is 
and the correlation histogram of emission from the cavity-centre system is shown in Fig.|5ja). 
The experiment assumes the use of a detector with a time bin size of 10 ps. The result shows a 
series of peaks separated by the clock period of the pulse. The suppression of the coincidences 
observed at zero delay signifies antibunching. More importantly, the suppression is observed 
during the period of a single excitation cycle for the short excitation pulse width T = 0.56 ps. 
The single photon probability per excitation trigger, estimated from the ratio of the number of 
single photon events to total number of pulses simulated is 0.99, while multi-photon probability 
is zero. With an effective lifetime of 70 ps, the system is capable of operating at an excitation 
rate of 10 GHz. However, to ensure all photon pulses are well-localized within their respective 
time bins, a bit-stream rate of 1 GHz is preferable. Finally, in Fig. |3b), we show an increase 
in multi -photon probability with increasing excitation pulse width. In agreement with the result 
from Fig.lltb), a considerably long simulation is needed to observe an appreciable multi-photon 
probability of ff{\)% for T ^ 10^ ps, the simulation was performed over a trajectory period of 
1.5 ms and resolution set to 10 ps. 

4. Conclusion 

We have studied the effect of a cavity on an NV^ defect centre in enhancing its spectral prop- 
erties for the purpose of single photon generation for quantum computing and communication. 
Assuming an atomic- vibronic NV model in single-mode cavity, we have shown that by cou- 
pling the centre strongly to a high-Q/V [Q ^ i^(10^ — lO""), V ^ X^] cavity that is resonant 
with the ZPL and with excitation using a sub-picosecond pump, the cavity-centre system is 
capable of issuing a photon of wavelength 638 nm with high spectral purity. We predict that the 
cavity-enhanced NV centre can have an effective lifetime of 70 ps and linewidth of 0.01 nm, in 
contrast with an unmodified centre's photoluminescence lifetime of 11.6ns and spectral width 
of 150 nm. Photons are emitted with near unit single photon probability of 0.99 while main- 
taining small multi-photon probability i^(10^^), thus making it a relatively efficient triggered 



photon source compared to a bare NV centre or an attenuated laser Finally, the device can 
potentially operate at a repetition rate of 1 GHz, considerably greater than demonstrated NV 
systems forQKD applications [32,33]. 
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